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Abstract 

We describe the infrared R-operation for subtraction of infrared 
divergencies in Feynman diagrams. 
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Massless particles can introduce infrared divergences (IRD) in Feynman 
graphs. In Minkowski space there are IRD in the S matrix for QED due to 
the emission of soft real photons, and also in graphs with virtual photons. 
These IRD cancel in the cross section [|[]. Furthermore, in QCD and also in 
QED in the approximation that quarks are massless, one has the situation 
that massless particles couple to themselves or to other massless particles, 
and this leads to further IRD in the S matrix in Minkowski space, the so- 
called collinear divergences. These cancel also if one averages over the color 
of the incoming particles and sums over the momenta of the initial states with 
soft collinear particles 0, or if one uses factorization methods. In this note 
we discuss IRD in Green's functions in Euclidean space; they are unrelated 
to the IRD in Minkowski space which occur in the S matrix. 

Four-dimensional quantum field theories with only dimensionless coupling 
constants contain in Euclidean space for generic external momenta no IRD. 
Thus in the proof of renormalizability of proper graphs in QED, all diver- 
gences which one encounters (in Euclidean space) are ultraviolet divergences 
(UVD). The Z factors contain thus only information about the small- distance 
behavior, and for this reason they can be used to construct running coupling 
constants. The same applies to QCD and massive quarks. Even if one sets 
the masses of all particles to zero, it remains true that for generic Euclidean 
external momenta the proper Green function in four dimensions are free from 
IRD[|. 

Spontaneously broken field theories have in general dimensionful coupling 
constants. For example, the 0(2) Goldstone model and the SU{2) Higgs 
model contain an interaction term Xva^ with Xv a dimensionful coupling 
constant. Nevertheless these four-dimensional Goldstone and Higgs models 
do not contain any IRD in the Euclidean proper graphs; this is due to the 
Goldstone theorem which states that proper selfenergies 11 (p) for Goldstone 
bosons vanish at = even when loops with massive scalars a contribute 
to the Goldstone boson selfenergy. 

However for such theories as Ay?^ + h(f^ theory in four dimensions with 
a superrenormalizable dimensionful coupling constant h, IRD do occur. A 

^In dimensional regularization one sets the following giuon selfenergy graph 
— Q to zero. This graph has an UVD but not an IRD, and its vanishing should 

be considered as the result of a computation, not as an independent rule. (One may for 
example replace -p- by + ^-^.^^^^^ , and one finds then that the sum of both inte- 
grals vanishes). Similarly — Q ^ J d^k/k^ vanishes, but now the IRD (evaluated 

at n > 4) cancels the UVD (evaluated at n < 4). Again this cancellation follows from the 
rules of dimensional regularization. 
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simple graph which shows this exphcitly is as follows 



k+q 




We have a massless boson in the larger loop with loop momentum k and a 
selfenergy insertion with loop momentum q due to a massive scalar. The two 
massless propagators /c~^ lead to a logarithmic IRD J d'^k/k'^ because now 
the proper selfenergy of the massive scalar Il{k) = J g2^jn2 (fc_g)2_|_^2 C^^g does 
not vanish for small k. 

Also at exceptional momenta IRD in Euclidean space can occur as the 
following example shows 




The two external lines in the middle carry momenta p and —p and one can 
also write them as a 93^ -insertion. The loop integral yields J d^k^ [k^qy 
with Q = Qi + Q2-, and contains clearly an IRD at A; = 0. 

Another area where IRD create problems is in the calculations of higher- 
loop (5 functions. Suppose we want to compute the divergences in the fol- 
lowing 4-point graph 




(3) 
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It simplifies the calculation a great deal if one sets ^3=^4 = 0. When we 
set some external momenta to zero, we shall call this operation "nullification 
of momenta". Nullification of some of the external momenta makes loop 
calculations a lot easier but it creates spurious IRD. (By spurious we mean 
that for generic momenta there are no such infrared divergences). To compute 
the UVD in this graph after having put Ps = P4 = 0, one must first remove 
the spurious IRD §]. 

The conclusion is that in four-dimensional Euclidean space Feynman 
graphs may contain IRD in the following cases 

(i) if there are superrenormalizable coupling constants such as \ip^ interac- 
tions in Z) = 4. 

(ii) at special external momenta. More precisely, when the sum of some ex- 
ternal momenta is equal to the sum of some (or none) of the other external 
momenta; in particular when some external momenta are nullified. (In the 
example with external momenta +p and — p, the sum of these two external 
momenta vanishes). 

When there are both massless and massive particles in the theory, the 
situation is much more complicated, and one must proceed by studying each 
theory separately. For example the interaction g^px^ in D = 6 with a massless 
field if and a massive field x leads to IRD in the selfenergy with two or more 
closed X loops,but if one first renormalizes and imposes the renormalization 
condition that the (p self-energy vanishes like fc^, the IRD disappears. This 
example shows that in general one should first renormalize and then 
study IRD. 

In higher spacetime dimensions the degree of IRD in general decreases 
because of the measure d^k, but in lower spacetime dimensions one en- 
counters more IRD. In particular, m D = 2 the nonlinear sigma models 
with action gijd^ip'^d^ip^ have dimensionless coupling constants (for example 
Qij = Siji^l+gip"^) has the dimensionless coupling constant g) but tadpoles and 
selfenergies are IRD. On the other hand, Xip^ theory in D = 6 dimensions is 
renormalizable but not superrenormalizable (because A is now dimensionless) 
and has no IRD.Q 

To explain these various results on IRD one can use some simple IR power 
counting. Consider a proper graph with nonvanishing external momenta pi 
in D dimensions. We shall assume that the external momenta are nonexcep- 

^As an example, consider in 6 dimensions a selfenergy graph with a massless scalar in 
the loop, and insert into this loop a string of M selfenergies with massive scalars in the 
loops. Then the propagators yield a factor (-p-)*^"*"^ and the measure yields / (fk, but 
now each massive renormalized selfenergy yields a factor J cPq (^qi_^_^n-i-^(\._qyi_^.jyj;i ~ ^ and 
there is indeed no infrared divergence. 
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tional by which we mean that there does not exist a relation 

M N 
j=l k=l 

for any M > and N > other than overall energy-momentum conservation. 
The propagators contain loop momenta and external momenta. Choosing a 
particular momentum flow through the diagram, there are "soft propagators" 
with only loop momenta and "hard propagators" with a combination of loop 
momenta and external momenta. (There are no propagators with only a 
combination of external momenta since the graph is proper). For vanishing 
loop momenta, the hard propagators do not become singular if we do not 
have exceptional external momenta. Hence, for nonexceptional values, the 
external momenta provide an infrared cut-off for the hard propagators. We 
can then determine the IRD which occur if one or more loop momenta tend 
to zero by shrinking all hard propagators to a point. The following example 
in A</?^ theory illustrates this procedure 




The contracted graphs are still proper when the original graph was proper. 
Because the hard propagators form a connected graph, there is only one 
contracted vertex for a given proper graph when the external momenta are 
nonexceptional. 

Just as one can count the degree of UVD of a proper graph by UV count- 
ing rules one can also develop IR counting rules. To perform IR counting, 
consider a contracted proper graph with N external lines, L loops, / internal 
propagators and vertices Vj with j lines, in addition to the contracted vertex. 
Let Ni be the number of soft lines at the contracted vertex (the subscript i 
stands for internal) . Because the contracted graph is still proper the number 
of soft lines connecting it to the rest of the graph is at least 2, hence Ni>2. 
Then the usual counting rule for the number of loops and the relation which 
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states that any internal (external) line ends at two (one) vertices, lead to 




N + 21 = ^3V, + {N + Ni) (6) 
j 

(In the example L = 4, = 6, V4 = 2, J = 6 and A^i = 4. Then L = 4 = 6-2 
and N +21 = 18 = 8 + 6+4). If there are no superrenormalizable couplings, a 
vertex Vj in D dimensions carries D — ^j{D — 2) momenta attached to it. (For 
example, in gauge theory in D = 4, the AAdA vertex carries one momentum 
and the AAAA vertex carries no momentum). When all loop momenta tend 
to zero, minus the overall degree of IRD of a contracted graph is given by 

uiR = DL-2I + J2Vj(^D-^jiD-2)^ 

= l{D-2)N,>D-2 (7) 

Therefore in D > 3 there are in general no overall IRD, but in D = 
2 all contracted graphs are logarithmically IRD. In particular the widely 
used WZWN models contain IRD. On the other hand, gauge theories in 4 
dimensions have no IRD, as we already discussed. 

What happens if only some of the loop momenta tend to zero, but others 
do not? If a number In ^ I of the internal momenta are kept hard, one can 
further contract the original diagram such that only the soft lines {I — In 
in number) remain. For example, one could make the loop momentum of 
the soft loop at the top in figure (|^) hard. If one were to make the two 
propagators in the loop on top hard, contraction of this new hard lines would 
yield a second contracted vertex, but contraction of one of the loops on the 
side would still leave only one contracted vertex. 

Let the new hard lines form a subgraph with Lh loops, Ih propagators 
and VjH vertices. The remaining degree of IRD is in this case 

uJiR (subcase) = ujir - Aujr 

AujR = DLh - 2Ih + (^D - ij(D - 2)^ (8) 

The subgraph can be worse IR divergent than the original graph because ujr 
(subcase) is equal to ujji of the original graph minus the degree of infrared 
divergence Aum of the subgraph. The original set of hard lines together 
with the new set of hard lines form a new (possibly disconnected) set of hard 
lines, and we can again apply the IR counting rules to the new contracted 
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graph. In the example in (|^) with 1 = 6 soft hnes we can make the two soft 
hnes on the right hard. Then the soft part subdivides as follows 

uJiR (subcase) = 3D - 8 + {-D + A) = 2D - 4 , . 

AuiR = D-4:+{D-2{D-2)) = 

The subcase remains indeed IR finite in D > 3. As another example, one 
may make the loop on top of figure hard. Than one finds that Auir = 
-D — 4 + {—D + 4) = —D + 4, and a;/i^( sub case) = 3D — 8, which is again IR 
for D > 3. However, an example of a subgraph which is divergent while the 
minimal graphis not is the following self energy in D=6: 



(10) 

The graph is overall IR finite (as we have proven generally) but contracting 
the two subloops, one finds an IRD proportional to J d^k/k^. We must thus 
learn how to subtract IR subdivergences. 

There exists a general scheme for subtracting IRD, which is the counter 
part of the so-called i?-scheme of BPHZ for subtracting UVD Q . The scheme 
which subtracts both UVD and IRD is called the R* scheme [^, and one can 
formally write 

R* = RuvRiR (11) 

According to this prescription, one first removes all IRD of given graph, and 
only afterwards subtracts all UVD. This is correct for most cases, but there 
are graphs where RuvRir is not equal to RirRuv- The question then arises 
which order is the correct one, and the answer is that the correct order is 
RirRuv^ and not (0). We already argued before that in the case of theo- 
ries with massive and massless particles one should first renormalize before 
extracting IRD. In most of our examples we shall follow the prescription 
in (^) because this is technically easier, but we shall also discuss a 5-loop 
graph where ([TT| ) is not correct. 

The UVD can be canceled by the usual UV counter terms, but the IRD 
are discarded by hand. For the calculation of (3 functions this is no problem 
because the IR divergences due to nullification of external momenta were any- 
how spurious, but for field theories such as h^^ 'm. D = A discarding genuine 
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IRD by hand seems a dubious procedure. One would prefer to have also an 
IR renormalization procedure similar to the UV renormalization procedure, 
but it does not seem to exist. Speculations have been made that the sum of 
infrared divergences vanishes in the two-dimensional WZWN model (when 
properly summed). 0] If this does not happen in this model, or in massless 
superrenormalizable theories, one would have to exclude such theories. 

In the R subtraction scheme of BPHZ, graphs are expanded into a Tay- 
lor series in the external momenta. We shall instead use dimensional reg- 
ularization to compute both the IRD and the UVD. Before going on, we 
make a comment on the relation J d^k/k'^ = in dimensional regulariza- 
tion. The reason this integral vanishes is that it contains both an UVD and 
an IRD, whose sum cancels. (One may separately compute the UVD from 
/ dt^k/^k"^ + rn?Y and the IRD from J ^^^^j^w^^ show that their sum 
cancels). If one were to use J d^k/k"^ = in higher-loop f3 function calcula- 
tions, one would drop some UVD, and hence one would make an error. In 
fact, one never encounters the need for setting J d'^k/k^ = in the computa- 
tion of the (3 function at the one-loop and two-loop level, but at higher loops 
care is required not to discard UVD. The claim is that the R* scheme does 
not loose UVD even though it sets tadpoles to zero according to the rules of 
dimensional regularization. 

Let us now explain the infrared subtraction procedure by some simple 
examples. Consider massless Ay?^ theory in D = 4. The proper graph 



contains no IRD at generic p, hence there is nothing to subtract. This is 
clear after contracting the graph 



Let us introduce an operator Rir which projects out the infrared finite part 
from a given graph. Then in this example we obtain 




(12) 




(13) 



RiR — e 



e 



IR-finite 



(14) 
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However, consider in D = A the graph 




(15) 



where the dot indicates that two external momenta have been nulhfied. There 
is then clearly a logarithmic divergence J df^k/k^. This IRD is due to the dou- 
ble propagator • — • — • and we can determine the divergence it produces by 
inserting it into the simplest graph where it yields an IRD. So we work in two 
steps: we determine first the IR divergent part associated with • • • and 
then we use it to subtract the IRD from the original graph. We indicate the 
IRD induced by the double propagator • • • ^■^ ( | ) ' ^^^^ ^ 1 ) 

should make the simplest graph with an • • • insertion finite as far as 
IR divergences go 




= IR-finite 



(16) 

Note that we obtain the remainder of the graph after extracting the IRD by 
deleting the subgraph which yields the IRD, not by contracting it. Analyti- 
cally the meaning of this graphic equation is as follows 

We can also formulate the subtraction procedure by replacing ^ in the 
original graph by a term with 5'^{k) 

^-^ + c,V(A:)/.%e = Z^-4 (18) 

The factor /i^ is needed in dimensional regularization to keep the dimension 
of the last term the same as that of fc"*^. Given the rules of dimensional 
regularization, one can compute Ci. 

We now return to the original graph in ([15|) . As before, the infrared 
finite part of the original graph is projected out by the operator Rir^ and 
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pictorially one has the following relation 

RiR ) = + il)m " IR-finite (19) 

Note that one always begins with the original graph on the right-hand side, 
and then one adds terms which subtract the IRD of the various subgraphs 
of the original graph. By convention one always writes plus signs on the 
right-hand side. 

The UV subtraction scheme can be formulated in the same way. Consider 
for example the graph in D = 4 




(20) 

It has a logarithmic UVD due to the loop in the middle, and also two loga- 
rithmic UV divergences due to the two 2-loop subgraphs on the left and on 
the right. (These two subgraphs are overlapping but that does not modify 
the subtraction rules). Let us introduce an operator Ruy which projects out 
the UV finite part of a graph. One obtains then 

Again by convention we always use plus signs for the terms to be subtracted. 

After the IRD have been subtracted, one may subtract the UVD. One 
does this for each term on the right-hand side of the Rjr equation separately. 
Consider for example (p^Qj). The subtraction of UVD proceeds as follows 

(I), o 

(22) 



uv 



IR 



We never subtract IRD or UVD from counter terms, so there is no "nesting" 
of the subtraction procedure. 

The UVD are easily located. Only the subgraph UVD. Hence 



= UV-finite 

uv 



(23) 

Ruv -O =-0 + (0)^^-^^ = UV-finite 
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The tadpole graph in the first hne vanishes according to the rules of dimen- 
sional regularization. (We set it to zero even though it contains an infrared 
and an ultraviolet divergence. The ultraviolet divergence in this diagram is 
accounted for by the whole R* procedure). Note that again we always begin 
with the original graph on the right-hand side, and then subtract divergences 
by adding counter terms. (Again by convention we use plus signs for these 
subtractions). 

Combining the IR and UV subtraction procedure, we obtain for the graph 



in (|T5D 



R* 



-e- 



uv 



-e- 



uv 



+ 



(1). o^(I)„(0 



(24) 



finite 



uv 



We can now determine the overall UV counter term which makes the graph 
finite after all sub divergences have been removed. This is the UV counter 

term one needs for the /3 function. It is given by ( — © — )uv and can be 

computed by evaluating the r.h.s. of the following equation 

^1 =-i^e- + (0; 



uv 



uv 



•+ 



(25) 



Taking the pole parts (PP), we can also write 



uv 



+ finite parts 



■e- 



uv. 



-e- 



pp 



(iL -o-(lL(0 







uv 



uv 



IT 



■+ 



(26) 



J pp 



Each graph on the r.h.s. is computed with dimensional regularization, in- 
cluding the original graph (which is of course the most difficult to compute). 
The subtraction terms denoted by {)uv and ()//? are polynomials in i, so 
one must compute also some graphs on the right-hand side to order e, etc. 
More precisely when there is a higher order pole ^ due to {)uv and {)ir, one 
must compute the corresponding graph to order e*^~^ 
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To illustrate the analogies and differences of Ruv a-nd Rir consider the 
following example 



+ (0)uv -Q- 



(27) 

The hatched part of the graph denotes any nonsingular subgraph. As this 
example shows, to remove UVD one shrinks the divergent proper subgraphs 

to a point, but one deletes the IR divergent subgraphs. We compute 

as the ^ pole part of this one-loop graph, and ( | ) is computed from 



R 



IR 



IR 



finite (noi) (28) 



In both cases one gets j poles (if one uses dimensional regularization) . 
Recall that the graphical identity 



R 



IR 



IR 



(29) 



corresponds to the following analytical expression 
1 1 



R 



IR 



{k - Qf 



/a (J, + _J_ (30) 



The IRD occurs at A; = and is taken care of by the ^5'^{k) insertion. One 
may find products of such IR factors 5'^{k) / e with difi^erent loop momenta, but 
never with the same loop momentum. As an example consider the following 
Feynman graph in ip'^ + ip^ theory 




dVqd^p 



{k - pYp'^k'' {Q - kf {Q - qy{p - qy 



(31) 
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In this example there is an IRD a,t k — p, and an overall IRD sd, k — p — 0. 
To subtract these IRD we replace some propagators by delta function in D 
dimensions 



Rjn{F) = F + / ^ f,^C,S^{k - p)) ^^j^ 



(Q-g)2 (p-g)2 



1 



iQ-k)HQ-gf {p-qY 



(32) 



Because there are now two D-dimensional Dirac functions one needs the 
factor /i^*^ with e = D — 4, as usual in dimensional regularization, to make 
the dimensions come out correctly. 

We can write this pictorially as follows 



R 



IR 





+ 



IR 




+ 



(33) 



-Q-^ = IR-finite 



IR 



We recall that we must first compute Ci and C2, and then by substitution 
we should find that the r.h.s. is IR-finite. 

The computation of Ci was discussed before, it follows from — . 

The computation of C2 follows from requiring that the following expression 
be IR-finite 



R 



IR 



4>)=4>*(i),. 




+2 



(34) 



• • • 



IR-finite 



* / IR 

The factor 2 is due to the two overall IRD, one upstairs and one downstairs. 
Note that (• • ) = f but f «^f^ ) _ = 7 + ^- the two lower prop- 



eran i 



IR 



agators would have been massive, we would have needed a factor one). Note 
also that • • • — Having determined Ci and C2 we can determine 

RlR . 

The IR counter terms are always polynomials in - and 6'^{kj). However, as 
our next example shows, one sometimes needs derivatives of 6'^{k). Consider 
the following massless graph 
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/ 



(35) 



P-Q 



There is now a quadratic IRD proportional to J d^p/p^. We want to 
extract it from the original graph by adding a term involving S^{p) but in 
order that dimensions match, we need a Dalembertian acting on 5^(p) 

d^p 



R 



IR 



p%p - QY 



d^p 



+ l^'CiDpS''{p) 



(36) 



Partially integrating the operator Dp 
the integrand 



{p - QY 

TTTTir- we obtain for the last term in 



d^pS^{p)Dp 



(p - Qy 



1 



(p - Qy 



(37) 



p=0 



To fix Ci we consider the simplest graph with this divergence; this is un- 
fortunately the graph itself but we nevertheless proceed (we could make the 
original graph more complicated to lift this degeneracy). Hence we fix Ci 
from 



R 



IR 





+ 



(38) 



The two slashes in the last part of the equation indicate the action of Dp 

2e(5~^, so we find the equation 



on l/(p — Qy. However, □p(p — Q)' 
* II * = — 2e • • • and 

1 1 



p=0 



« / + Ci2e—- = IR- finite 

p'{p-Q? Q^ 



(39) 



One would expect that Ci is proportional to ^, so the original graph happens 
to be IR- finite (due to the peculiar properties of dimensional regularization) . 

To evaluate Ci we therefore need another graph. One could take a massive 
propagator with {p — QY in which case 

'8-2D){p-Qy -2Dm^ 



0,[(p-QY + m']-' 



[{p - QY + m 



213 



(40) 
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is no longer proportional to e. One can then determine Ci. 

We now discuss a subtlety having to do with the order in which one 
applies Ruv and Rir. The combined operation is denoted hj R*. A priori 
one might expect that RirRuv is equal to RuvRm, but there exist counter 
examples at the five- loop level 0. Consider the following 5- loop graph 




Operating with RjjvRir gives the incorrect result 
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+ 



+ 



.. (jSl)-[jgl.(0)..-@L 

(5>)..^-(^l 

(i).[-«{-^-<0L-^-(2>)..^ 

{^*(^)J 



2 € 



+ 8 



+ 



+ 



(^),.[^*(^)J 

(^l[o*(^)J 



(42) 

The IRD arise when the following loop momenta vanish: p, pql, and pqls. 
When pt vanish (or pi), one finds an IRD of the term J dpdtp~'^t~'^ , but this 
was already considered in the IRD with J dpp~^ so we do not count it sepa- 
rately. We repeatedly used that tadpoles vanish, no matter how many loops 
they contain. We also used (^)-(|39|), but because there are two propagators 
on which Dp can act, we get also cross terms where each propagator carries 
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one derivative. Operating with RmRuv gives the correct result 



•(O)„(0)..[^.(^).— ] 

■(oi(o)..[^-(^).. — ] 
•(*).. 

(43) 

As another example, we consider a 5-loop graph which is needed to compute 
the /3 function at 5-loops. The original diagram gives a vertex correction for 
the A</7^ coupling, namely a graph with 4 external lines, but we nullify all 4 
lines. Since vacuum graphs vanish in dimensional regularization we add two 
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new external lines carrying a new external momentum Q. 





Q 




(44) 

The reason we let momenta Q flow in and out the graph at these particular 
points is that this allows to compute the original graph easily. Indeed, to 
compute the graph itself (which is always needed), one may first compute 
the subgraphs 



2e 



2+2 e 



(45) 

Then one evaluates / d^q{Q — qf) ^^(?^) ^ We now evaluate R* on 
this diagram. Since there are no ^ terms, there is no term with \I\p5^{p) 
and thus no ambiguity whether one should choose RjjyRjji or Rj^Ruy. 
We choose the former. We first record the result for acting with Rj^ on 
the graph, and then record the results due to acting with Ruv on each of 
the terms in the result for Rir. We write the results such that each col- 
umn in the result for Rjjv corresponds to one term in the result for Rir. 



R 



IR 




+ 2 




-OO- 



The subgraph ^ / ) (which corresponds to the IR divergent inte- 

gral J d'^kd'^p{k — p)~'^k p^^ for small k and p) is written as c6^{k)6^{p) 
whereas the subgraph ^ j contains only one delta function 5^ (fc-|- 

q) (since there is only one external momentum, one can only use one S^{p)). 
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At the end no IRD are left ( -QO- is IR finite). 

Now we perform Ruy on each of these terms. We write the result of each 
Ruv operation as a column. 



(i) 



{ ^ 



\ 



<DO- 



(46) 



Only the term at the bottom of the first column is needed for the (3 function. 
Since we need all \ poles, we must calculate each graph to some power in 
e. In particular, since the IRD in the second column is proportional to ^ we 
need the graphs in this column to order e°, but since the IRD in front of the 
third column is ^ + ^, we need the graphs in this column to order e, and we 
need the terms proportional to l,e and in the graphs of the last column 
since the IRD in front of the last column contains a leading term The 
three 1-loop, 2-loop and 3-loop tadpoles graphs in the first three columns in 
the third row all vanish according to the rules of dimensional regularization. 

All these calculations were done using dimensional regularization, not 
dimensional reduction (the latter is inconsistent at higher loops [§). For 
theories with 75 and t^upa the approach followed is to first compute graphs 
without e tensors (substituting for 75 the product ^^upal'^l^l^l")) and only 
after the calculation is finished one contracts with four-dimensional e tensors. 
In fact, to compute RuvF of a divergent graph with several e tensors, one 
can use the fact that RuvF is UV finite, and write the product of two e 
tensors in terms of D-dimensional 5 functions. The error one makes is of 
order e, so vanishes as e — > 0. These D-dimensional Kronecker delta function 
one can then insert inside the expression RuvF to obtain a scalar. If one has 
a single e^yp^-, one can multiply by kfk2k2kleai3-ys and work out the product 
^afi-yS^^'^^" in terms of D-dimensional delta functions. Effectively this means 
contracting the open indices in RuvF with momenta k^ to obtain a Lorentz 
scalar. So, in the end one never computes with open indices. One obtains the 
correct answer for the 3-loop chiral anomaly This approach works only 
for multiplicatively renormalizable quantities, and not diagram-by-diagram. 
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The reason is that for multiphcatively renormahzable models 

RuvF = ZF (47) 

Then the error in using D-dimensional contractions is of order D — A. 

Another example where the correct subtraction of IRD is crucial for de- 
termining the UVD is the massless WZWN model in D = 2 dimensions 
The simplest IR subtraction corresponds to 

RiR -Q- = -Q- + 2 ( J)^^ • • = IR-finite (48) 
A more complicated example is 

We write the numerator in terms of the momenta which appear in the 
propagators, qp{k — 2q)„ = qp{k — q)„ — qpqa, and obtain then graphically 

(50) 

The slashes denote momenta, and reduce the IRD. We obtain then 

(51) 

Having subtracted the IRD, one may then proceed to compute the UVD. 
These are supposed to cancel in the D = 2 WZWN model at L > 2 loops, 
but one clearly needs to be careful with first subtracting the correct amount 
of IRD. 
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qp(k - 2q)„ 
{k — qYq'^{k — pY 



d^qcPp 



(49) 



If there are no momenta in the numerator, one obtains 



-©-=-0- *Kl)„-C)-*^(0),/ — ' = 

(52) 

If one of the propagators is massive, we obtain (denoting the massive prop- 
agator by a solid line) 



(53) 

We close with two further examples. First consider in D = 4 the following 
2-loop graph in massless Xip'^ 

(54) 

Suppose wc want to compute the 2-loop contribution to the Z factor of the 
vertex. We first nullify two external momenta because this simplifies the 
calculation 




Next we subtract the IRD 



R 



IR 



= iR finite 



(56) 
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Finally we subtract the UVD 



""(-©-•(ll^)--a-(o)„-j 

•(^),-(il-^-(l).(oL 



uv ^ ♦ ^IR 

(57) 

The tadpole graph vanishes, and the contribution to the Z factor follows then 
from 

(^l-[^-(ll--(ll(o)J,,„, 

(58) 

For fun we give one last example |10]. We consider the following 3-loop 
graph in massless A<y9^ in D = 4 




a 

UVD — Q— = 

b 



\ ; -Q- = 

f 



a 



b 



The infrared subtraction yields 



(59) 



R. 



(60) 

Note the appearance of a disconnected graph. Next ultraviolet subtraction 
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yields 



(o) O 

(O) -£l 



(O) (O) _Q 



(I), 




■(^);((0i.- 



The UVD are UVDi:^^) = UVD{^2] 



are IRD (73 



-1 

(47r)2e 



(61) 

with € — n — A, while the IRD 



and IRD{^4) is of the order 4- The overall UVD of 



this graph is Pq = "3^3(4^)6 (where G denotes the original graph), and R*G 
is UV and IR finite. 

We end with final conclusions and comments 

(i) Graphs in which all masses have been set to zero allow one to compute 
[3 functions in a much simpler way than keeping masses, but one introduces 
spurious IRD which one should subtract. We have explained the rules for 
subtracting IRD from Feynman graphs and for determining the final UVD 
which are relevant for f3 functions. 

(ii) All UV counter terms and IR counter terms are polynomials in -. Tad- 
poles can be set to zero; even though they may contain UVD; in the total 
result for the UVD as given by the R* scheme no UVD are lost. 

(iii) One can write the counter terms as Feynman graphs with some prop- 
agators replaced by 5^{k) or \3k5^{k). For example, in D = 2 one has 
RiR J ^ = / d'^k(j2 + 7^^(^))- In this sense the IR counter terms are local 
in p space. In x-spacc they would be nonlocal, because infrared divergences 
deal with the large x behaviour. 

(iv) To remove UVD one shrinks subgraphs, but to remove IRD one deletes 
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subgraphs. The final UVD is the one one needs for j3 functions. For lower 
loops RjjvRiR is equal to RirRuv, but for higher loops the order matters, 
and the correct order is R* = RjrRuv- The order only matters if a graph 
contains a factor 

(v) in the original BPHZ approach |lT|, one starts with RuvF = (1 — 



t^)nj7e$(l — t^)F where H are all proper subgraphs of the Feynman dia- 
gram F which are superficially divergent (power-counting divergent), while 
yields the overall divergence after all sub divergences have been subtracted. 
Furthermore, if H D H' one should write 1 — to the left of (1 — ), but 
if H and H' are disjoint or overlapping, it does not matter in which order 
they appear. There exist refinements which show that one only needs subsets 
of subgraphs which are nonoverlapping ("forests"). One can then prove the 
forest formula RuvF = (1 — t^) J2i^He<t,i{—t^)F where the 0, are a forest 
which includes the empty set. In all examples above we have evaluated this 
forest formula both for UVD and IRD. 
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